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^ ! Abstract: We derive a compact expression for the Borel sum of a QCD amplitude in 

terms of the inverse Mellin transform of the corresponding Borel function. The result al- 
lows us to investigate the momentum-plane analyticity properties of the Borel-summed 



D , Green functions in perturbative QCD. An interesting connection between the asymp- 



totic behaviour of the Borel transform and the Landau singularities in the momentum 

plane is established. We consider for illustration the polarization function of massless 

rS I quarks and the resummation of one-loop renormalon chains in the large-/3o limit, but 



p^ ', our conclusions have a more general validity. 
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1. Introduction 

It is well-known that the moinentuin-plane analyticity properties imposed by causality 
and unitarity are not automatically satisfied by the correlation functions calculated in 
perturbative QCD. This is actually a more general feature of approximate solutions 
of physical equations. For instance, unexpected analytic properties were discovered 
for the perturbative solutions in potential theory [[![]. In QCD, the problem is more 
severe due to confinement. Indeed, as shown in [Q, in a confined theory the poles 
and branch points of the true Green functions are generated by the physical hadron 
states in the unitarity relation, and no singularities related to the underlying quark and 
gluon degrees of freedom should appear. On the other hand, in perturbation theory the 
branch points of the amplitudes correspond to the production of free quarks and gluons. 
Other, more complicated singularities make their appearance if one goes beyond naive 
perturbation theory. For instance, renormalization-group improved expansions have 
unphysical space-like singularities (Landau poles or cuts), while the introduction of 
non-perturbative contributions such as vacuum condensates in the Operator Product 
Expansion (OPE) produces poles and branch cuts at the origin ||^. Both types of 
singularities are absent from the exact QCD amplitudes. 

Recently, the contribution of large orders of perturbation theory related to the so- 
called renormalons [^ have been investigated by several authors, using in particular 
the method of Borel summation |§]-[Q (see [^ for a review and further references). 
In general, these resummations introduce additional momentum-plane singularities in 
QCD Green functions, reflecting the presence of the ultraviolet (UV) and infrared (IR) 
renormalons in their Borel transforms. But despite of the large number of works de- 
voted to Borel summation in perturbative QCD, the momentum-plane analyticity of 



the Borel sum did not receive particular attention so far. The reason presumably lies in 
the fact that perturbation theory by itself is not complete and must be supplemented 
by non-perturbative terms to build up the exact theory. It is therefore not surprising 
that theoretical predictions based on perturbation theory are plagued by bad analyt- 
icity properties absent from the true physical amplitudes. Nevertheless, given that 
analyticity is the key property allowing the continuation of theoretical predictions from 
Euclidian to Minkowskian kinematics relevant to hadronic observables, we believe an 
investigation of the momentum-plane analyticity of the Borel sum in perturbative QCD 
is of both theoretical and practical interest. 



The importance of the problem was emphasized some time ago by Khuri p8[, who 
noted the mathematical possibility of a Borel sum with analytic and asymptotic prop- 
erties in the momentum plane different from those of the individual terms in the per- 
turbative expansion. (Khuri also quotes unpublished results by J. P. Eckmann and 
T. Spencer illustrating this possibility.) The deep connection between momentum- 
plane analyticity and Borel summability was revealed by 't Hooft ^, who showed that 
the presence of resonances and multi-particle thresholds on the time-like axis is in 
conflict with the mathematical conditions required for the existence of the Borel sum. 
More recently, the relation between the asymptotic behaviour of the Borel transform 
and the momentum-plane singularities of the QCD amplitudes was discussed in [|l^] . 
Renormalization-group invariance and the Landau pole in the running coupling asiQ"^) 
are the essential ingredients of this analysis. Attempts to remove this unphysical sin- 
gularity in the context of dispersion theory with an IR-regular running coupling were 
made in [^. This procedure raised much interest, since it generates power corrections 
of a type not present in the OPE ||3D|, ^ (see, however, the discussion in [0). Further 
clarifications of the connection between the Landau singularity of the running coupling 
and the IR renormalons of the Borel transform are desirable in this context. 

In the present paper, we discuss some aspects of the momentum-plane analyticity 
of Borel-summed Green functions in perturbative QCD. To this end, we consider the 
polarization function for massless quarks and calculate its Borel sum at complex values 
of the momentum transfer. The result we derive, written in a compact form in terms 
of a "distribution function" (i.e., an inverse Mellin transform of the Borel function), 
allows us to investigate the singularities generated by the process of Borel summation. 
Since we consider the massless case, branch points at s = will be present in the 
time-like region, which of course do not reproduce the true hadronic thresholds. Our 
main concern, however, is the appearance of unphysical (Landau) singularities in the 
space-like region. Our discussion is kept at a rather formal level, although reference to 
the specific case of a single renormalon chain (corresponding to the so-called large-/?o 
limit |TB|, |2^) is made to illustrate various results and assumptions. In Section 0, we 
start with a brief review of the integral representations of the Borel function in terms of 
its inverse Mellin transforms. Two different types of inverse Mellin transforms (called 



distribution functions) were introduced in 15 and 18 . We indicate the mathemati- 



cal conditions under which these functions exist, the relations among them, and their 
connection with renormalons. An important finding is that the inverse MeUin trans- 
form of the Borel function associated with the polarization amplitude is, in general, a 
piece-wise analytic function, with two pieces defined inside and outside a circle in the 
complex plane, each of them being analytic except for a cut on the real axis. In Sec- 
tion 1^, we derive a compact expression for the polarization function at complex values 
of the momentum variable in terms of the distribution function defined in [J15l . The 



same technique is applied for obtaining an expression for the Borel sum of Minkowskian 
quantities, using as input the standard result for the Borel transform of these quantities 
obtained by analytic continuation. In Section ^, we discuss the analyticity properties 
of the polarization amplitude in the momentum plane. These properties depend on the 
prescription used for treating the IR renormalon singularities of the Borel transform. 
We adopt a principal-value prescription, which leads to amplitudes without space-like 
singularities outside the "Landau region" (the region between and the location of 
Landau singularities on the space-like s axis in the renormalization-group improved 
perturbative expansions of correlators in QCD) and is therefore consistent with the 
analyticity requirements imposed by causality and unitarity The analytic continua- 
tion to low momenta reveals explicitly the connection between the Landau singularities 
and the piece-wise character of the inverse Mellin transform, which in turn is deter- 
mined by the asymptotic behaviour of the Borel function. Our conclusions are given 
in Section |^. The paper has an Appendix, in which we present a method for deriving 
relations between the various distribution functions based on Parseval's theorem for 
Mellin transforms IBl . 



2. Borel transforms and distribution functions 

Consider the current-current correlation function 11 (g^) defined by 

-AttH I d'x e'"?-" (0| T {V^{x), ¥"{{])} |0) = {q^q"" - g^^'q^) n{q^) , (2.1) 

where V^ = q'j^q is the conserved vector current for massless quarks. We define the 
corresponding Adler function 

D{s) = -s^, s = q\ (2.2) 

as 

which is UV finite. The function n(s) can be obtained from D{s) by logarithmic 
integration, yielding 

s 

n(s) = - /d In(-s') D{s') + const. , (2.3) 

where the integration is along a contour in the complex s' plane ending at the point 
s and not encountering the singularities of the integrand. General causality and uni- 
tarity properties of QCD Green functions imply that n(s) and D{s) are real analytic 



functions in the complex s plane (i.e., n(s*) = n*(s) and D{s*) = D*{s)), cut along 
the positive real axis from the threshold for hadron production at s = Am^ to infin- 
ity. The analyticity properties of the exact Adler function D{s) are summarized in the 
dispersion relation 

AM ^£/d/ ■■""<-' y . (2.4) 

4m| 

where the spectral function Imn(s + ie) is non-negative. The function n(s) satisfies a 
once-subtracted dispersion relation similar to ( |2.4| ). 
Writing the one-loop running coupling constant as 

An 
"^^-^) = /3oln(-s/A2)' ^'-'^ 

where Po = ^^c — f'^/ is the first coefficient of the /3 function and A the QCD scale 
parameter, one finds that the first terms in the renormalizat ion-group improved per- 
turbative expansions of the functions D{s) and n(s) are given by 
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U{s) = k - \n{-s/A') )-^-^ + . . . , (2.6) 

where k is a. constant. These expressions are derived using the renormalization-group 
equations in the perturbative region |s| 3> A^. Their analytic continuation to low 
energy contains unphysical singularities on the space-like axis, which are absent from 
the exact amplitudes. Specifically, D(s) has a Landau pole at s = — A^, and n(s) has 
a Landau cut along the interval — A^ < s < 0. In |^ it was proposed to restore the 
correct analyticity properties by a redefinition of the coupling ( |2.5D using a dispersion 
relation. In the present paper, we take a different perspective and analyse the impact 
of large orders of perturbation theory on the analyticity properties of the correlation 
functions. 

The structure of high-order contributions to the perturbative expansion of the Adler 
function was investigated by many authors. The expansion coefficients are known to 
exhibit a factorial growth, indicating that the perturbation series has zero radius of 
convergence 0, |^. If one formally applies the Borel method to sum the series, this 
growth is refiected in singularities of the Borel transform on the real axis. Let us 
write the formal Laplace integral expressing D{s) in terms of its Borel transform in the 
large-/3o limit as [^ 



oo . ^,s 
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l + ^yd«(^) Bniu), (2.7) 



where C is a scheme-dependent constant, which takes the value C = —5/3 in the MS 
scheme and C = in the V scheme, and Ay = e~^ h? is a scheme-independent combina- 
tion of this constant with the QCD scale parameter. The corresponding representation 
for His) is 



oo 



— S \ 1 /■ , —s 



nW = *-'"UJ-^A/^"UJ B„(„). (2.8) 

We have used the renormalization-group invariance to isolate the dependence on the 

momentum variable s, so that the Borel transforms Bj:){u) and Bji{u) are scale and 

scheme independent and also do not depend on s. Note that the definition (|2.2| ) implies 

the relation 

B^{u) = -Bd{u). (2.9) 

u 

The Borel transforms Bd{u) and Byi{u) have singularities on the real axis in the 
complex u plane. Some of them, the so-called IR renormalons, are situated along the 
integration contour in ( p.7|) and (12.^ ).-^ The precise nature of these singularities is only 
known for some specific cases. For instance, the first IR and UV singularities are branch 
points located at m = 2 and u = —1, respectively, with a universal nature determined by 
renormalization-group coefficients [^ ^. Information on the other singularities is only 
available from calculations performed in the large-/?o limit, which predict the correct 
locations of the singularities but not their nature (in general, the large- /3o approximation 
yields pole singularities rather than branch points). Specifically, the expression for the 
Borel transform of the Adler function is [^ § 

g,,.,) , _lji- g (-1)^* .^. (2.10) 

3(2 - u) t, [k^ - (1 - uf-f 

The presence of IR renormalons renders the perturbation series not Borel summable 
and the Laplace integrals ill defined. In order to evaluate them, a prescription for 
avoiding the singularities must be specified. A question we address in this paper is 
whether additional information can be employed to prefer one choice of a prescription 
over another. We shall argue that momentum-plane analjd;icity might be a helpful 
guiding principle in this respect. 

We also consider Minkowskian quantities such as the finite-energy moments of the 
spectral function of n(s) defined as 



,fc- 



- f dss''lmU{s + ie) = -j^ (fdss''U{s), (2.11) 



4ml \s\=so 



where Sq > 4m^ is an arbitrary cutoff, and we have applied the Cauchy relation for the 
true, physical polarization function n(s) to pass from the integral along the cut to a 



-'^The function i?n(u) has, in addition, a pole at u = related to the short-distance behaviour of 
the current correlator. This pole is removed by standard UV renornialization. 



contour integral along a circle in the complex s plane. Combinations of such moments 
enter the prediction for physical observables such as the hadronic decay rate of the 
r-lepton [^, and some tests of QCD related to the dependence under variation of Sq 
are discussed in [^, ^. The Borel representation for Minkowskian quantities such as 
the spectral moments Mk can be derived starting from the contour integral written in 
the second relation in (|2.11| ) and inserting the expression ( p.8|) for n(s). The result is 



M. 
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Bm^{u) 



du 






Bm^ (u) , So > A 



2 
V ; 



k+ 1 



smvrM 



Bd(u) 



(2.12) 



(2.13) 



k + 1 — U TXU 

Next we define the "distribution functions" coresponding to the above Borel trans- 
forms. We consider first the Borel transform Bd(u) of the Adler function and assume 
that it is analytic in the strip u\ < Reu < U2, where ui = —1 and U2 = 2 are the 
positions of the first UV and IR renormalons, respectively. Let us further assume that 
the following L^ condition holds: 



uo+ioo 

— - / dn|5i5(n)p < oo, 

UQ—ioo 



(2.14) 



where Uq G [mi, U2\- Then the inverse Mellin transform of Bj:){u) can be defined as p4 

<1uBd{u)t''~\ (2.15) 



Wd[T] 



1 

2Td 



The condition ( |2.14| ) is satisfied by the Borel transform calculated in the large- /3o limit, 
given in ( |2.1CI| ). In this case, the function wnir) was first introduced in [^, where its 
physical interpretation as the distribution of the internal gluon virtualities in Feynman 



diagrams was pointed out. Relation (|2.15|) can be inverted to give 



Bd{u) = jdTWD{T)T~\ 




and the following completeness condition holds p4 

UQ+ioo 
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2^ 



&a I-Bd(^) 



dr 



t^-^''wd{t)\\ 



UQ—lOO 







(2.16) 



(2.17) 



The integral relation ( |2.16| ) defines the function B£,{u) in a strip parallel to the imagi- 
nary axis with wi < Reu < M2. 

The distribution function idnlr) can be calculated from (|2.15|) by closing the inte- 
gration contour along a semi-circle at infinity in the u plane and applying the theorem 



of residues for the singularities of Bij{u) located inside the integration domain. For 
r < 1 the contribution from the semi-circle at infinity vanishes if the contour is closed in 
the right half of the u plane, while for r > 1 the contour must be closed in the left half 
plane. One thus obtains different expressions for the distribution function depending 
on whether r < 1 or r > 1, which we shall denote by w\j (r) and id\j (t), respectively. 
By Cauchy's theorem, these functions are given by the contribution of the residues of 
the IR and UV renormalons, respectively. For instance, wj^ (r) has the expression 



where the contours C+ and C_ are lines slightly above and below the positive real axis, 
respectively. A similar expression applies for the function w\j (r) , where the integration 
lines are now parallel to the negative real axis. In the large- /3o limit, these functions 
are given by ||15| 



WD\r) = y |r (^^ - Inr ) + (1 + r) |L2(-r) + Inr ln(l + r) 



— U+nr+- + -nr- 
3l V4 2 J T 



w'^^\t) = ;;U + lnr+ (; + ;inr)- + 



L2(-r-i)-lnrln(l + r-i)]| 



(2.19) 



where L2{x) = — Jq — ln(l — t) is the dilogarithm. The above expressions are analytic 
in the complex r plane, with no singularities other than branch cuts along the negative 
real axis. Together they define a function u}£)(r) that is piece-wise analytic in the cut 
r plane, with different functional expressions for |r| < 1 and |r| > 1. 

The fact that the function wd{t) is piece-wise analytic will have interesting con- 
sequences in our analysis below. One might ask whether these analytic properties are 
valid in general, i.e., beyond the large-/3o limit. It is important in this context that we 
consider the analytic continuations of the function w\) (r) defined in (|2.18 ) and of the 



function tD|) (r) defined in a similar way, and not the original definition ( |2.15| ). The 
equivalence between the two results is valid only for real r. By writing r" = exp(M In r), 
and noticing that the residues of the renormalon singularities are real, we obtain from 



( |2.18|) real values for positive r. However, a branch point at r = and a cut at negative 



r may appear due to the logarithm. In general, the singularities of the Borel transform 
are expected to be branch points (rather than poles) located on the real u axis. The 
inverse Mellin transforms of functions with branch points have a similar structure as 
those of functions with pole singularities |^: they are piece- wise analytic functions, 
composed of two pieces analytic in the domains |r| < 1 and |r| > 1, with a cut along 
the negative real axis. Therefore, this seems to be the most general analytic structure 
for the function wd{t), assuming that the asymptotic behaviour of the Borel transform 
is in accord with the norm condition (p.l4|) . 



The application of the above procedure for the polarization function itself requires 
some care due to the additional pole at u = of the Borel transform Bii[u). An 



expression consistent with the definitions (2^) and (2^) is obtained if we define the 
inverse Melhn transform wu{t) in terms of an integral along the line Keu = Uq with 
tio > 0. Then, as proven in the Appendix, the result is 



T 

■u^nlT) = — dx wj:,{x) . 



(2.20) 



In the large- /So limit, the explicit expression for this function is [^ 
... 16 ' 



X 



wnVT) 



1 — In X H — (5 — 3 In x) + 



1 + x 



X 



L2(-x) +lnxln(l + x) ), (2.21) 



where a: = r if r < 1, and x = 1/r if r > 1. It follows that w^i {t 



(»i 



w 



«)i 



1/r). 



Consider now the case of Minkowskian quantities such as the spectral moments 
Mfc. In this case, the presence of the factor sinvru in the Borel transforms in ( p.l3| ) 
affects their asymptotic behaviour for large u and invalidates the L^ condition ( |2.14|) . 
Therefore, in general an inverse Mellin transform defined as in ( |2.15| ) does not exist for 
Minkowskian quantities. This observation suggests to extract the sine factor from the 
Borel transform and consider the inverse Mellin transform of the remaining expression. 
Actually, such a definition can be considered not only for Minkowskian quantities but 
for any Borel transform Bx{u) fl^. We define a function 



hx{u) = 
and assume that it obeys the L^ condition 



Bxiu) 



siianu 



(2.22) 



no+ioo wo+ioo 
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Bxiu) 



smvr-u 



< CXD 



(2.23) 



where uq ^ 0. In the case of Minkowskian quantities, the sine function in the denomi- 
nator of ( ^.23| ) compensates the factor of sin ttu appearing in the Borel transform. (In 
the case of Euclidean quantities the sine factor brings an additional improvement in the 
convergence of the integrals.) Provided the condition ( ^.23| ) is satisfied, we can then 
define the inverse Mellin transform Wx{t) as 



Wx(t) 



1 
2^ 



UQ+ioo 



1 r , Bxiu) 



dubx(u)T^-' = — I dM ''^'^^-^-1 



2TTi 



■ r 



SmTTM 



(2.24) 



This definition was introduced in a more physical context in |T8[. Relation 
be inverted to give 

bx(u)= IdrWxHT 



241) can 



(2.25) 



which together with (|2.22|) defines the Borel function Bx{u) along a fine Keu = Uq > 0, 
where Uq can be taken arbitrarily small. Of course, if the Borel transform Bx{u) itself 
satisfies an L^ condition like in ( |2.14| ), then also the normal inverse Mellin transform 



wx(t) defined as in ( |2.15| ) exists. In this case, one can prove the relation [^ 

feM = iL5£M. (2.26) 

TT J X + T 



which we derive in the Appendix using Parseval's theorem. The representation ( p.26| ) 
shows that, unlike the distribution functions W£)(t) and w-air), which are piece- wise 
analytic in the r plane, the modified inverse Mellin transforms Wd{t) and Wii{t) are 
real analytic functions in the entire complex r plane cut along the real negative axis. 
(The explicit expression for the function Woij) can be found in the Appendix of ||IB[.) 
On the other hand, for the Minkowskian quantities the distribution functions Wx{t) 



are piece-wise analytic, with two pieces defined for |r| < 1 and |r| > 1 |18[ . 



3. Borel sums in the complex momentum plane 

Borel-summed expressions for the polarization function at both space-like and time-like 



momenta in terms of the distribution functions wd{t) and Wd{t) were derived in JTS 



and [|I^, respectively. In what follows, we will generalize the techniques used in these 
works to the calculation of the polarization function in the complex momentum plane. 
For definiteness, we adopt as a regularization prescription for the ill-defined integral 
( p.7|) the so-called "principal value" , which amounts to taking one half of the sum of 
the integrals along two parallel lines slightly above and below the real axes, denoted 
by C± (the same integration lines appeared previously in ( p.l8|) ). We thus define 



D{s) 
where 



d+{s) + d- 



(3.1] 




d±{s) = ^ I du{-^\ Bd{u), (3.2) 



and s is taken to lie outside the Landau region, \s\ > Ay, so that the integrals are 
convergent. The prescription ( |3.1| ) is a generalization of the Cauchy principal value for 
simple poles, which was adopted in |3^ as the only choice giving a real value for the 



Borel-summed amplitude when the coupling constant is real. In the renormalization- 
group improved expansion, this means that D{s) is real along the space-like s axis 
outside the Landau region. This is consistent with the dispersion relation (p^.4D, which 



follows from general causality and unitarity requirements. As we shall see, other pre- 
scriptions violate these requirements. 

Our aim is to express the integrals (|3.2| ) in terms of the distribution functions defined 
in the previous section. To this end, we must pass from the integrals along the contours 



C± to integrals along a line parallel to the imaginary axis, where the representation 
( p.l6| ) is valid. This can be achieved by rotating the integration contour from the real 
to the imaginary axis, provided the contour at infinity does not contribute. In each 
individual case, one must establish that this condition is satisfied. Let us consider 
first a point in the upper half of the momentum plane, for which s = \s\e^'^ and 
—s = \s\ e*^'^"'^^ with a phase < < tt. Taking u = TZe^^ on a large semi-circle of 
radius TZ, the relevant exponential appearing in the integrals (|3.2|) is 



exp < —TZ 



In — 2" cos 9 + {n — (p) sin 6 

At/ 



(3.3) 



For |s| > Ay, the exponential is negligible at large TZ for cos9 > and sin 6' > 0, i.e., 
for the first quadrant of the complex u plane. Assuming that Bd{u) increases slower 
than any exponential [|I^, the integration contour defining (i+(s) can be rotated to 
the positive imaginary axis, where the representation ( ^.l(j| ) is valid. This leads to the 
double integral 



«oo oo 

d+{s) = — du /drw£)(r)exp 
Pn J J 



-u lln^ + i(0-7r) 



(3.4) 



The order of integrations over r and u can be interchanged, since for negative {(f) — n) 
the integral over u is convergent. Performing this integral yields 



d. 



wJ'^ 



WD{r 



H~rs/Al 



(3.5) 



Consider now the evaluation of the function d-{s) given by the integral along the 
contour C_ below the real axis. Naively, one might think to rotate the integration 
contour to the negative imaginary axis without crossing any singularities. However, 
this rotation is not allowed, because along the corresponding quarter-circle sin 6' < 0, 
and the exponent ( |3.3| ) does not vanish at infinity for (p < tt. The way out is to perform 
again a rotation to the positive imaginary u axis, for which the contribution of the 
circle at infinity vanishes. But in this rotation the contour crosses the positive real 
axis, and hence we must pick up the contributions of the IR renormalon singularities 
located along this line. This is easily achieved by comparing the expression (|2.18| ) for 
the function tD}^ (r) with the definition ( p.2|) of the functions d±{s). It follows that 
d-{s) can be expressed in terms of d+{s) as 



(i_(s) = d^{s) 




'^)w^n^\'Al/s) 



(3.6) 



The relations (|3.5| ) and ( |3.6D completely specify the function D{s) in the upper half of 
the momentum plane. Using the same method, the Adler function can be calculated in 
the lower half plane, which corresponds to taking tt < (p < 27r. In this case, the integral 



10 



along C_ can be calculated by rotating the contour up to the negative imaginary u 
axis, while for the integration along C+ one must pass across the real axis. Combining 
the results, we obtain the following expression for the Adler function in the complex 
momentum plane: 



DW(,) = 1 + ^ dr ^^Sl^ T ^ -^ w'i}\-MM . (3,7) 




where the superscript "±" in parenthesis refers to the sign of Ims. The corresponding 
expression for the polarization function n(s) can be obtained by inserting the above 
result into the definition (|2.3| ). This gives 




oo 






nW(.) = A:-ln -^ --^/dru.,,(r)lnln -|iU^ -^ 4<^(-Ay^) 



(3.8) 
with wu{t) as defined in (|2.2CI|) . Note that the last terms in (|3l7| ) and (|3.8| ) involve the 



analytic continuation of the Mellin transforms ^^(r) and wu{t) from the real positive 
axis, where they have been calculated, to arbitrary complex arguments. 

Before discussing in the next section the analyticity properties of these results, 
we apply the same technique to the principal-value Borel summation of Minkowskian 
quantities. Let us generically denote a Minkowskian quantity by R{sq), with sq some 
fixed scale, and write its Laplace integral with the principal-value prescription as 

Riso) =Ro + \ [r+iso) + r_(so)] , (3.9) 

where Rq is a constant, and 



r±{so) = — — / dubR{u) I exp 
2tPn Jc+ 




U In —rr + ilXU 



exp 




u In —rr — iTCU 



^3.10) 
We have expressed the Borel transform Bji[u) in terms of 6_r(m) as in (p.22|) and com- 
bined the two exponentials arising from the factor sinvrw with the exponential in the 
Laplace integral. We recall that sq > Ay. The procedure described in detail for com- 
plex s can be now applied in a straightforward way. For C+ (C_) one can rotate the 
integral of the first (second) term in (|3.1CI| ) towards the positive (negative) imaginary 
axis, where the integral representation ( |2.25| ) of b^lu) in terms of Wr{t) can be used. 



For the remaining terms, i.e., the second (first) term in the integral along C+ (C_), 
the simple rotation towards the imaginary axis cannot be performed, since the cor- 
responding terms ±i7ru in the exponential blow up along the quadrants of the circle. 
As explained before, we must first cross the real axis in the u plane and then make 
a rotation towards the nearest imaginary semi-axis. But in crossing the real axis we 
encounter the IR renormalons, whose contribution must be added using an analog of 



11 



(3.6). In the present case, the relation is 



/ dubR{u) <exp 
/ dubniu) lexp 



—win 



—win 



So 

' So ' 



± inu 



± inu 



HO 



2mT±W}i^>iT±), (3.11) 



where t± = —Ay/sQ^ie. Using this relation, we find after a straightforward calculation 



R{s„) 



Po 



dr 



W.ir) 



In^TSo/AD + n^ 



+ 



2/3o So 



KO 



«)/ 



W^A"^(r+) + l^r(^- 



(3.12) 



Notice that the last two terms involve the values of the function Wj^ (r) for real 
negative values of the argument, where this function has a cut. In the two terms the 
cut is approached from opposite directions, and the sum gives twice the real part of the 



function W 



R ' 



r) at r 



'Ay/ So- It is convenient to rewrite the above result using an 



integration by parts, yielding 



1 
R{so) = Ro + -^ drlji^r) arctan 
Po -l 



IT 



ln(rso/A^) 



^I/SO 



+ ^Re / drlR{T-te), (3.13) 

On "' 



-Al/so 



where 



Ur) = - 



d 
d7 



rWjiir 



(3.14) 



The integral representations ( |3.12| ) and ( |3.13| ) were derived previously in |]T8| by means 
of a different technique of treating the integrals (|3.10|) . Specifically, in order to avoid 
the crossing of the real axis in the rotation of the contours C± towards the imaginary 
axis, these authors use formally the integral representation ( p. 251) of bji{u) outside its 
range of validity. The double integral thus obtained is evaluated by performing first 
a suitable rotation of the integration line in the r plane, and then a rotation of the 
integration contour in the u plane. The formal application of this procedure leads to 
an expression identical to ( |3.12| ). 

For our further discussion, it will be useful to have an explicit expression for the 
functions Im^. (r) entering the integral representation ( p.l3| ) in the particular case of the 
Borel sum of the spectral moments M^. These functions are given by 



hi At) 



k + 1 



7r 



dxx woir/x) 



k + l 



■ T 



TC 



dz z 



-k-l 



Wd{Z 



r > 0. 



(3.15) 



We stress that from the derivation of the result (|3.13|) it follows that for negative 
values of r the functions Im^{t) must be obtained from the analytic continuation of the 

^In this reference, results are given for the functions Wkir) — ttt iM^i''')- 
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expressions derived from for r real and positive. Note that for r < these functions 
are not the same as the functions 



I'm^t) 



k+ 1 



7r 



-t)'' Uzz-^-^wd{-z) 



r < 



(3.16) 



one would obtain using the analytic continuation of the function wj:){t/x) under the 



integral in ( |3.15|) . This observation will become important in the next section. 



4. Analyticity properties of the Borel-summed functions 



Let us now investigate in more detail the momentum-plane analyticity properties of 
the Adler function and of the polarization function calculated in the previous section. 
We start with the Borel sum of the function D{s). From (|3.7|) , it is apparent that 
its s dependence is isolated in two terms having a different origin. The first term can 
be regarded as an average over gluon virtualities of the one-loop expressions in (|2.6| ) 
1^]. It is obtained using the integral representation of the Borel function in terms 



of the inverse Mellin transform along the imaginary axis in the u plane. The second 
term is generated by the residues of the IR renormalons when crossing the real axis to 
pass from one of the contours C± to the other. As we will see, the fact that this term 
contains the piece w\) (r) of the distribution function has important consequences on 
the analyticity properties. 

As shown in ( ^.T] ), we obtained two different expressions for the Adler function valid 
in the upper and lower half of the momentum plane. Using general theorems on the 
analyticity of functions represented by integrals 



^«)^ 



I, and the fact that w\) (r) is 



holomorphic for complex values of r, it follows that the functions D^^\s) are holomor- 
phic for complex values of s outside the real axis, in the upper and lower half planes, 
respectively. We will now show that they represent a unique function analytic in the 
cut momentum plane. To this end, we apply the Schwarz reflection principle. The func- 
tions D^^\s) are related to each other by the reality condition [D'^^\s*)]* = D^~\s), 
since w}j (r) is a real analytic function, which enters the expressions for D^^\s) and 
D^~\s) with imaginary coefficients of opposite sign. Assume now that s is approach- 
ing the real axis from the upper half plane. Using expression ( |3.7|) together with the 
relation 



Im 



ln[— (x + ie)] 
we obtain on the time-like axis 



9{x) 



vr 



vr 



luiD{s + ^e) = -5- 

Po 



dr 



WD(r 



In'(x) 



A' 



n 



TC^ 



6{x + l) 



In 



2Vs/A2 



y. 



+ 77^ 



+ '-^Rew^<\-Al/s) 



(4.1) 



s>0, (4.2) 



while on the space-like axis 




lmD{s + ^e) = ^ wo{~K/s) - w'^'i-A^s) 



i<)i 



s < 0. 



(4.3) 
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We recall that we have assumed that s lies outside the Landau region, i.e., |s| > Ay. 
In this case, the argument of the functions ^/^(r) and w)^ (r) appearing in ( [4.3| ) is 
less than one. But for such values the two functions coincide, and the two terms in 
( [4.3| ) compensate each other. Therefore, the imaginary part of the Borel-summed Adler 
function D{s) vanishes for space-like momenta outside the Landau region: 

lmD{s + ie) = 0, s< -A^ . (4.4) 

The same result is obtained when approaching the real negative axis from below. There- 
fore, in the deep Euclidian region —s> Ay both expressions in ( ^^71) lead to a real value 
given by 

D{s) = l + -RefdT .""^7,,. , .<-A^ (4.5) 

/5o ^ ln(-rs/Af) 

By the Schwarz reflexion principle, the two expressions D*-^^(s) in ( ^.71 ) define a single 
function D{s) real analytic in the momentum plane, with no cut for s < —Ay. 

Before proceeding, it is worth emphasizing that the vanishing of the imaginary 
part on the space-like s axis outside the Landau region is not a generic feature of the 
Borel sum for the Adler function, but is specific to the principal-value prescription 
adopted above. It does not happen if one takes a different prescription for avoiding the 
singularities of the Laplace integral. For instance, taking the integral along the contour 
C+ as a definition of the Borel integral we would obtain instead of ( |3r7| ) the expressions 



/?o J H-rs/ Al 







Di-\s) = 1 + 1 /dr-^?^ + ^ ( -^ 1 w^^\-A-ys) , (4.6) 




(3o{^' H-Ts/ Al 

where the subscript "+" indicates the regularization prescription. It is obvious that 
these expressions do not satisfy the condition of real analyticity D*{s) = D{s*). In 
particular, using the relation ( |4.1j ) we see that the two expressions in ( [4.6| ) have non- 
zero imaginary parts along the whole space-like axis s < 0. (These imaginary parts 
coincide, while the real parts exhibit a discontinuity across the real axis.) On the other 
hand, for real analytic functions the discontinuity is due to the imaginary part, since 
in this case D{s + ie) — D{s ~ ie) = 2ilm.D{s + ie). By unitarity, the discontinuity of 
the exact function D{s) is non-zero only above the threshold for hadron production, 
and below this threshold the exact function must be real.^ The expressions in (|4.6| ) are 
in conflict with the above requirements, and in our opinion this fact is an argument in 
favour of the principal- value prescription. 

The explicit expression for D{s) derived in ( |3.7| ) for |s| > Ay can be analytically 

continued inside the circle \s\ < Ay. It is interesting to make this continuation in 

^These properties, which follow from the general principles of field theory, are explicitely incorpo- 



rated in the dispersion relation (2^), due to the fact that the spectral function Iinn(s) is real. 
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order to explore what happens, after Borel summation, with the Landau pole present 
in the fixed-order perturbative expression ( |2.6| ). We find that for complex values of s 
outside the real axis the resulting expression is holomorphic. The discontinuities along 
the real axis can be analysed using ( [4 .21 ) and ( [4 .31 ). On the time-like axis, we find a 
rather complicated behaviour besides the expected unitarity branch point at s = 
produced by the logarithm. Indeed, the spectral function at s < Ay involves the values 
of the function Wj^ (r) extended analytically in the range |r| > 1, where there may be 
additional singularities (see, e.g., the explicit expressions (|2.19|) valid in the large-/?o 
limit). Singularities are also present along the space-like region —Ay < s < 0. Indeed, 
in this case the argument of the functions wd{t) and w]-, (r) in ([4.3|) is Ay/|s| > 1. 



Therefore, the last two terms in ( [4.3| ) do not compensate each other anymore, and hence 
the analytic continuation of the expression for D{s) derived above has a cut along the 
low-energy Euclidian region with the discontinuity 



\mD{s + ie) = — ^ 




w^^\-Kl/s) - w^^\-Kl/s)\ , -A^ < . < . (4.7) 

The appearance of this cut is explicitely related to the fact that the inverse Mellin 
transform of the Borel function is a piece-wise analytic function. Note that even if IR 
renormalons were absent and the theory was Borel summable, the Landau cut would 
persist. Indeed, in this case d+{s) = d_{s), and from ( p.6|) it follows that the function 
Wjy (r) vanishes identically. But in ([4.7|) we are then left with the term wjy (t), which 
receives contributions from the UV renormalons. 

The presence of the Landau cut ( [4. 7] ) is intimately related with the asymptotic be- 
haviour of the Borel transform. It is straightforward to show that if the Borel transform 
vanishes sufficiently fast at infinity in the complex u plane - more precisely, if not only 
iSoiu) but also the product Bd{u) sinvrw are L^ integrable functions in the sense of 
( p.l4|) - then the Landau cut is absent. Indeed, in this case a reasoning similar to that 
presented in the discussion containing equations ( p.22| )-( p.26| ) leads to the conclusion 
that the inverse Mellin transform u)d(t) satisfies a dispersion representation and is 
therefore analytic in the complex r plane cut along the negative axis. Then the two 
pieces w^ (r) and Wj^ (r) coincide, and the imaginary part ([4.7D vanishes identically. 



However, the strong decrease of the Borel transform ensuring the absence of the Lan- 
dau cut is not observed in the large- /3o limit, and it is thus rather improbable that it 
should happen in the physical case. 

As a side remark, let us indicate what happens if the effective coupling constant is 
modified in the IR domain so as not to contain a Landau pole. To study this case we 
make the replacement 



M ( 47r 



-^ ^ rz. = exp , (4.8) 



where the coupling as{—s) is defined in the V scheme. Then the discontinuity (|4.3| ) 
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across the Landau cut may be written as 



Im D{s + le) = — TL\wD{rL)-w'i^\rL)] , s<0. (4.9) 



The couphng as{—s) may be defined by a dispersion integral [^ or simply by sub- 
tracting the pole at s = —Ay by hand. As discussed in [Q, this redefinition amounts 



to adding terms exponentially small in the coupling, which do not modify the pertur- 
bative expansion. The modified coupling is real and positive along the whole space-like 
region s < 0, implying that tl defined in ( [4. 81 ) is a positive number less than unity. As 
a consequence, the imaginary part of D{s) vanishes along the space-like axis, since for 
Tl < 1 the two terms in ([4.9| ) compensate each other. Therefore, unphysical singular- 
ities do not appear in the Borel-summed expansion in powers of a regular coupling, if 
the principal-value prescription is adopted for treating the IR renormalons.^ 

For completeness, we also present results for the polarization function n(s). The 
Borel-summed expression (p.8|), obtained for |s| > Ay, can be analytically continued to 
the whole complex plane. We note that 11 (s) is holomorphic for complex values of s and 
satisfies the reality condition n(s*) = n*(s). On the real axis, this function can have 
singularities manifested as discontinuities of the imaginary part. A straightforward 
calculation gives 



OO / 

Imn(s + ie) = Tc + — drt(}£)(r) arctan 



■K 



ln(rs/A 



2 



-|- — Re / drM}/)(r), s>0, 

^0 A 2 / 

-Al/s 



vr f A^ 



Imn(s + ^e) = ^ I -^ 1 €jn{-Al/s) - idi^\-Al/s)] , s<0. (4.10) 

Po V S / L 

Again, the second expression vanishes outside the Landau region, where the function 
wu{t) coincides with w^ (r). However, inside the interval —Ay < s < the function 
n(s) has a non-zero imaginary part. 

The final point we will discuss concerns the implications of the above results for 
the Borel summation of Minkowskian quantities such as the spectral moments M^. In 
( |3.13|) we derived the principal-value Borel sum for these quantities, which we shall 
denote by M^™'^'. On the other hand, starting from the explicit expression ( |3.8| ) for 
the polarization function n(s) in the complex plane, we could calculate its moments 



also by performing the contour integral in the second relation in ( |2.11|) . Let us denote 



by M"''^^'^ the result of this procedure. We are now going to show that M^™''' = M™'^'*', 

^This argument cannot be considered a proof, since the very meaning of the Borel summation of 
a modified expansion with the Landau pole removed is not transparent. Indeed, the redefinition of 
the coupling itself may be seen as a partial summation of some higher-order effects, which must be 
separated from the terms taken into account in the Borel sum in order to avoid double counting. 
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which is equivalent to the statement that contour integration and Borel summation 
commute with each other. In order to compare the two resummed expressions for the 
moments, we pass from the integral along the circle |s| = Sq to an integral along the 
real axis. Applying the Cauchy relation, and taking into account the singularities of 
the Borel sum of n(s) inside the circle, we write in the most general way 

circle V^O K + 1 



M', 



27risg+^ 



ds 5*^0(5) 



\s\=ri 



+ 



k + 1 

IT Sq 



-V 



ds / Im U(s + ie) + / ds s'' Im n(s + ie 



-so 



(4.11) 



The point s = must be treated separately since the function Wn (—Ay/s) entering 
the expression for n(s) in ( |3.8| ) is singular at this point. (In the large-/5o limit, we have 
wfj (r) ~ r In r for r — > oo.) From this behaviour, it follows that the circle of radius t] 
gives a non-vanishing contribution only for the first two moments. For k > 2, the only 
singularities are due to the discontinuity of the imaginary part of n(s) given above. 
As we discussed, this imaginary part vanishes for s < —Ay. By a straightforward 



calculation, using (^4.10|) and isolating the contribution of the region —Ay < s < Ay, 



r circle 



we write ([4.1 1|) as Mj: 

^_^o k + 1 



27ri s 



k+1 



M'p. + 6k, where 
dss'^Uis) 



\s\=V 



+ 



k + 1 
Po^o 



a2 



I dss'=A(-A^/s)+ fdss''ReA{-Al/s) 



-A2 



with 



A( 



w 



(» 



(r) 



w 



«), 



r 



(4.12) 



(4.13) 



The quantity M'^ has an expression that is formally identical with the Borel moments 



^Borei jj^ ( p.i3| )^ but in which the functions Im,,{t) in ( |3.15| ) are replaced with the 



functions I'm^'^) i^^ ( |3.16| ) obtained using a different prescription for negative r. In 
other words, we have M^ 





^fc^"'"' - ^'k with 



TT 



^'k = lT^^ 

Po 



dr 



'A/. T 



It 



u^ 



le 



(4.14) 



Hr/so 



6'^, and we believe this result 

JVf Borel ^ 



An explicit calculation in the large-/3o limit shows that 5k 

is of general validity. From this observation, it follows that indeed M™^^^ 

which is in agreement with the findings of [|I^]. Note that the arguments given above 

show that the equivalence of the two resummation procedures results from a subtle 

compensation of the contribution of the singularities of n(s) along the Landau region 

with an additional term arising from the analytic continuation of a piece-wise analytic 

function. 
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5. Conclusions 

We have derived a compact expression for the Borel-summed Adler function and the 
polarization amphtude in the complex momentum plane in terms of the inverse Mellin 
transform of the corresponding Borel functions. At present, these inverse Mellin trans- 
forms (or distribution functions) are known only in the large- /5o limit. However, the 
expressions derived in the present work are useful for investigating the analyticity prop- 
erties in the complex momentum plane, and therefore some conclusions can be drawn 
even if the distribution functions are not known exactly. 

Our main results are the expressions for the Borel sums of the functions D{s) and 
n(s) given in ( p.7| ) and (p^), respectively. They have been derived using a generalized 
principal- value prescription, which leads to real values of the functions along the space- 
like s axis outside the Landau region, in accordance with general analyticity require- 
ments imposed by causality and unitarity. The expressions derived in the momentum 
plane outside the Landau region admit an explicit analytical continuation to the region 
|s| < Ay, revealing Landau singularities more complicated than those present in fixed- 
order perturbation theory. We have shown that the discontinuity across the Landau 
cut is connected with the piece-wise character of the inverse Mellin transform, which is 
in turn correlated with the asymptotic behaviour of the Borel transform. We have also 
presented a condition on the asymptotic behaviour of the Borel transform that would 
ensure the vanishing of the Landau cut. This condition is, however, not satisfied in the 
large- /3o limit, and it will most likely not be satisfied for real QCD. Therefore, we do 
expect the Landau cut to be a general feature of Borel-summed perturbation theory in 
QCD. 

With the explicit expressions derived in the complex momentum plane, we have 
checked that the spectral moments defined by contour integrals of the Borel sum of the 
correlator n(s) coincide with the "standard" Borel moments considered previously in 
[0]-pO|. This equality is consistent with the complicated singularity structure of the 
polarization function discovered in this work, providing an indirect check of our results. 

We have shown that the unphysical discontinuity across the Landau cut would per- 
sist in the fictitious case of a Borel-summable perturbation series. Thus, the Landau 
singularities in the momentum plane and the IR renormalons are not directly related 
to each other. A similar conclusion was formulated in the context of a specific model 
in [^ . The unphysical cut is a manifestation of the incompleteness of any perturba- 
tive approximation to a physical hadronic quantity. On the other hand, our results 
suggest that, if the running coupling is regularized by a dispersion relation, then the 
analyticity properties of fixed-order perturbation theory are preserved by Borel sum- 
mation. Indeed, in this case the discontinuity across the Landau cut vanishes even if IR 
renormalons are present. This result holds under conservative assumptions about the 
properties of the inverse Mellin transform and thus seems to have a general validity. 
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Appendix 

Here we indicate a simple method for calculating the distribution functions of various 
observables using the relations between their Borel transforms. To this end, we use 



Parseval's theorem for the Mellin transform |3J], which can be stated as follows: if 
J-'{u) and G{u) satisfy the L^ conditions 



uo+ioo l~uo+ioo 

-^ / du\J^{u)\'<oo, -^ / du\g{u)\'<oo, (1) 

cm J zm J 



where < mq < 1, then 

Uo+ioo oo 

— du T{u) g{l-u)= fdx fix) g{x) , (2) 

zm J J 

UQ—ioo 

where f{x) and g{x) are the inverse Mellin transforms 

UQ—ioc 1— UQ—ioo 



As a first application of this result, we use it to prove relation ( |2.26|) . We start with 



relation ( p.24|) and notice that the function Wx{t) can be written as 

uo+ioo 

W^x(r) = — / duJ'iu)gil-u), (4) 

uo— ioo 

where 

T(u) = B,(u), 5(u) = -£—^. (5) 

It then follows that f{x) = idx{x), and inserting the definition (|^) of the function Q{u) 
into (^ and making the change of variables s = 1 — t we find 

uo+oo -, _ 

9{x) = —- ds — . (6) 

zm J smvrs 

UQ—ioo 
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This integral can be calculated closing the contour in the complex s plane and applying 
the theorem of residues. For Re (t/s) > we close the contour in the left half plane 
Res < 0, while for Re (r/s) < we close it in the right half plane Res > 0. In each 
case, we pick up the contributions of the relevant poles situated at s = ±n with the 
residues 

ti): (I)" . (7) 

TIT \X/ 

The contributions of the semi-circles at infinity vanish due to the exponentials in the 
sine function, even for complex values oir/x = \t/x\ exp{±iilj) with \il)\ < vr. Therefore, 
we obtain 

,w^^|(-ir(i)"^^. (8) 

Inserting the expressions for the functions f{x) and g{x) into the right-hand side of 
(H), and using the definition of Wx{t), we obtain 

TT J X + T 



where r can take arbitrary values, except for real negatives. 
As a second example, we derive relation (|2.20|) , starting from 



uo+ioo g 

ffi„M = J- / d„M!i),.-., (10) 

2m J u 

UQ—ioo 

where we have used (|2.9|) . The above result can be rewritten as in (§), with the 
identification 

J'{u)=Bniu), g{u) = j^^. (11) 

A straightforward calculation gives the inverse Mellin transforms 

f{x)=wn{x), g(x) = l9(l-^y (12) 

Inserting these functions into the right-hand side of Parseval's theorem (|^), we obtain 
relation (|]20D. 
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